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The electromagnetic form factors provide important hints for the internal structure of the nucleon 
and continue to be of major interest for experimentalists. For an intermediate range of momen- 
tum transfers the form factors can be calculated on the lattice. However, reliability of the results 
is limited by systematic errors due to the required extrapolation to physical quark masses. Chi- 
ral effective field theories predict a rather strong quark mass dependence in a range which was 
yet unaccessible for lattice simulations. We give an update on recent results from the QCDSF 
collaboration using gauge configurations with Nf = 2, non-perturbatively 0(a)-improved Wilson 
fermions at very small quark masses down to 340 MeV pion mass, where we start to probe the 
relevant quark mass region. 
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1. Introduction 



In recent years the phenomenological interest in the electromagnetic form factors of the nu- 
cleon has revived. This was triggered by the Jefferson Lab polarisation experiments ^ ^ measur- 
ing the ratio of the proton electric to magnetic form factors, jj.^G < f\Q 2 )/Gm\Q 2 ). From these 
measurements an unexpected decrease of this ratio has been found, which means that the proton's 
electric form factor falls off faster than the magnetic form factor. 

Many theoretical calculations have been done to investigate possible interpretations of a de- 
crease of this ratio (see, e.g., [||] for an overview). Lattice techniques allow the calculation of the 
form factors from first principles. Such calculations do not only yield phenomenologically interest- 
ing quantities such as magnetic and electric charge radii and magnetic moments. These techniques 
also allow, e.g., the investigation of the Q 2 dependence of the nucleon electromagnetic form factors, 
which can be compared with experimental results and also helps in understanding the asymptotic 
behaviour of these form factors. 

In practice, the calculation of these form factors on the lattice remains a challenge. In recent 
years progress has been made to improve control on systematic errors which are related to the fact 
that the calculations are performed on finite volumes, at finite lattice spacings and at quark masses 
which are still relatively large. For making reliable predictions at physical quark masses it turned 
out that numerical results at smaller quark masses are crucial. With recent advances in computing 
power available to lattice QCD calculations and the speed-up of algorithms for simulating dynam- 
ical fermions it is now possible to reach much smaller quark masses with pseudoscalar masses in 
the range of 300 MeV. 



2. Calculation details 



In this talk we present results obtained on configurations with two mass degenerate flavours of 
non-perturbatively 0(a) -improved Wilson fermions. We choose Wilson glue for the gauge action. 
To scale the lattice results for different /3 and K" sea we use the Sommer parameter r (/3, K" sea )/a and 
the conversion factor ro = 0.467 fm to translate our results into physical units. 

The form factors are obtained from the standard decomposition of the nucleon electromagnetic 
matrix elements 



{p',s'\J^\p,s) = u(p',s' 



7^i(<2 2 ) + /a 



2M-, 



-F 2 (Q 2 



N 



u(p,s), 



(2.1) 



where p (s) and p' (s') denote initial and final momenta (spins), q = p' — p the momentum transfer 
(with Q = —q 2 ) and the nucleon mass. By calculating the matrix elements on the l.h.s. and 
the nucleon mass we obtain the Dirac form factor F\ (Q 2 ) and Pauli form factor F2(Q 2 ). 

The nucleon matrix elements are extracted from ratios of three- and two-point functions: 



R(t,r,p, P )= x 



C 2 (T,p')C 2 (t,p')C 2 (t-Z,p) 
C 2 (T,p)C 2 (t,p)C 2 (t-T,p>) 



1/2 



(2.2) 



Here t denotes the location of the sink. Assuming the source being located at time slice 0, we 
expect an plateau for <C T <C t (source and sink are separated by a distance of ~ 1.1 fm). For 
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Figure 1: The iso-vector Pauli form factor radius at (j3, K" sea ) = (5.25,0.13575) (left plot) and the ratio 
of the iso-vector d- and w-fiavour Dirac form factors (right) as a function of Q 2 . In the right plot we show 
results for similar mps ~ 400 MeV but different lattice spacings. 

further details see [||]. We use three different polarisations r unpo i = ^(1 + 74), Ti = ^(1 + 74) /7571, 
^2 = |(1 + 74) '7s7> as well as three different sink momenta po = (0,0,0), p\ = (p,0,0), P2 = 
(0,p,0) (where p = 2n/Ls). 17 different choices for the momentum transfer q have been used. 
Due to statistical fluctuations the operand of the square root in Eq. ( |L2| ) may become negative. 
Results for which this happens are discarded from the consecutive analysis steps. 
We use the local vector current, which needs to be renormalised and improved: 

Vp =Z Y (\+bvam q ) [qy^q + icvadx(q~Opxq)] . (2.3) 

The renormalisation coefficient Zy and the parameter by have been determined non-perturbatively 
M\. The improvement coefficient cy is only known perturbatively. However, since this coefficient 
is expected to be a small number and because the improvement term was found to be small in the 
quenched approximation [j|], we will ignore the improvement of the operator. 

In the following we will consider the iso-vector, iso-scalar and the proton form factors. The 
latter two might receive contributions from quark-line disconnected terms, which are notoriously 
hard to calculate on the lattice and will not be considered here. 



3. Parameterisation and Q 2 dependence 

First we will investigate the Q 2 dependence of the Dirac and Pauli form factors. Both lattice 
and experimental data can be reasonably well parametrised by a pole ansatz 

= (31) 

From naive dimensional counting one would expect the Dirac and Pauli form factors to scale 
differently, i.e. F\ oc Q 4 and F2 °< Q 6 , which corresponds to p = 2 and p = 3, respectively. 
Experimental data as well as theoretical calculations indicate deviations from this naive picture. 
For instance, the JLab results for ^J&F^f' ] (Q 2 ) /Ff' ] {Q 2 ) were found to be surprisingly flat. A 
perturbative QCD analysis of the Pauli and Dirac form factors predicts the ratio 

{Q 2 /log 2 ^I^Q 2 /A 2 )F 2 {Q 2 )/F Y (Q 2 ) (3.2) 
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(with j3 = 1 1 — 2A/f/3) to scale as a constant [Qj. In an investigation of the Q 2 dependence of the 
experimental nucleon form factor data using empirical parameterisations, Diehl and collaborators 
found indications for the form factor scaling to be flavour dependent [f7j|. By comparing the Dirac 
form factor results for proton and nucleon they found the flavour contribution F-j® (Q 2 ) to decrease 
faster with Q 2 than f/ u) (Q 2 ). 

From a first inspection of the lattice results one finds that the effects from changing p in 
Eq. < ]3 - 1[ ) are small with respect to the statistical errors. In Fig. (left plot) we show the results of a 



fit to f[ v) for one particular data set using p = 2 and 3. It is clear that it is difficult to obtain lattice 
data with high enough precision over a large enough range of Q 2 values to distinguish between 
a dipole or tripole behaviour. It may, however, be instructive to consider ratios of form factors 
in order to reveal significant deviations from the naive scaling hypothesis. In the right plot of 
Fig. [I] we show the ratio Fj /Fj and find that it does not scale as a constant. This is consistent 
with the observation by Diehl et al. We should however emphasize that we are ignoring possible 
contributions from disconnected terms. 

Based on the above observation we perform fits to the Dirac and Pauli form factors for each 



flavour separately using the ansatz Eq. ( p.l[ ). To eventually decide which p should be used, we 
perform these fits for various 1.5 < p < 4 and search for the "optimal" p for which # 2 /d.o.f. 
becomes minimal. We observe strong variations of the resulting "optimal" p for data sets which 
differ in (/3,fC sea ,K" va i)- However, we nevertheless see clear trends. In case of the w-flavour Dirac 
form factor F^ (Q 2 ) we find for most data sets p to be close to 2. On the other hand, for the other 
form factors F^iQ 2 ), F 2 (u) (<2 2 ) and F 2 (d) (g 2 ) we typically get p ~ 3. 

4. Form factor radii and magnetic moments 

From the fits to Eq. ( |3.1[ ) we obtain the form factors at zero momentum transfer F,-(0) and the 
di- or tripole masses M,-. Equivalently, we can use the same fit to obtain the form factor radii r,- and 
the magnetic moment fi. These quantities are defined as follows: 



Fi(Q 2 ) = F(0) 



l-lr 2 Q 2 + 0(Q 4 ) 



(4.1) 



H = Fi(0)+F 2 (0). (4.2) 

From the magnetic moment we can calculate the anomalous magnetic moment, e.g. fc^ = ju,W — 1. 
For comparison with phenomenological results we use the normalised anomalous magnetic mo- 
ment, e.g. K"M norm = k"M m^{m K ) /m^(mps), where m^(mps) refers to the nucleon mass calcu- 
lated on the lattice at the quark mass corresponding to the pseudoscalar mass raps and m^(m n ) the 
experimental value of the nucleon mass. 

In Fig. H we show results for and F^ (0) as a function of nip S which have been calculated 
for different values of the gauge coupling /3 and various sea quark masses. These results seem 
to lie on a universal curve, which indicates that discretisation errors are small. In the following 
we will consider them as negligible compared to the statistical errors. We furthermore observe 
that the iso-vector anomalous magnetic moments and di-/tripole masses show a linear quark mass 
dependence for a very large range of quark masses. However, using an ansatz linear in (roraps) 2 
to extrapolate our results to the physical point, we obtain values which differ significantly from 
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Figure 2: The iso-vector Dirac form factor dipole mass (left plot) and the iso-vector Pauli form factor at 
zero momentum transfer (right plot) as a function of nip S . 

those extracted from experiment. This is consistent with calculations based on a chiral effective 
field theory (ChEFT) that includes nucleons, pions and delta resonances as explicit degrees of 
freedom [g[ g[ These calculations predict for the iso-vector form factor radii and the iso-vector 
anomalous magnetic moment a strong quark mass dependence in the small quark mass region. This 
region is just starting to become accessible for simulations with dynamical Wilson quarks. 

In the left plot of Fig. |3] we compare the lattice results for the iso-vector Dirac form factor ra- 
dius with the following result from ChEFT, where we used the same phenomenological parameters 
as in 01: 

,2 1 



(4^)2 

2 



l+7gi + (l0gi + 2)log [^]| 



+ 



CA 



547l 2 F 2 



26 + 30 log 



+ 30 



X 
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PS 




For the iso-vector Pauli form factor, rlj , and the anomalous magnetic moment, K^ v \ we per- 
formed a combined fit of the lattice results to the following expressions from ChEFT: 
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Figure 3: Lattice results for the Dirac (left) and Pauli (right) radii. The solid lines show the ChEFT results 
given in Eqs. (4.3) and (4.4). Note that the left curve is not based on a fit. In the right-hand plot the dashed 
line is the result of a fit to the tripole masses linear in (romps) 2 . The star denotes the experimental value. 



Here we keep the chiral limit of the anomalous magnetic moment fo v )°, the iso-vector N-/S. coupling 
cy and the ChEFT parameters B C 2 and as free parameters. The result of this fit is displayed in 
the right plot of Fig. |3] for (Vj ) 2 and the left plot of Fig. || for together with the lattice data 
and phenomenological results at mps = m K . 

The lattice data for the Dirac radius do not seem to agree well with the ChEFT result. Since it 
is not clear up to which quark masses the ChEFT expression is valid, results at even smaller quark 
masses will be needed to actually clarify this issue. For both the Pauli radius and the anomalous 
magnetic moment the lattice and ChEFT results look consistent. It is somewhat surprising that this 
seems to hold also for rather heavy quarks. From our preliminary results at very low quark masses 
we see first indications for the Pauli radius to bend towards the phenomenogical value. 

Finally, in the right plot of Fig. ^, we show our results for the iso-scalar form factor. From 




0.5 1 
nip, [GeV] 




Figure 4: Lattice results for the iso-vector (left) and iso-scalar (right) anomalous magnetic moment. The 
solid lines in the left plot show the result from a fit to the ChEFT expression in Eq. (4.5). In the right-hand 
plot an ansatz linear in (romps) 2 has been used to fit the data. The star denotes the phenomenological value. 
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ChEFT a linear quark mass dependence is expected, which is fully consistent with our lattice 
calculations. 

5. Conclusions 

We have presented the current status of the calculation of the electromagnetic form factors of 
the nucleon by the QCDSF collaboration. At the currently achieved level of statistical errors, still 
large uncertainties remain for the parameterisation of the form factor results. However, qualitative 
agreement with the experimental data has been found, e.g. the flavour dependence of the Dirac 
form factor F\ . As new configurations at very small quark masses are starting to become available, 
we are improving our control on the extrapolation of the lattice results towards the chiral limit. We 
have found first indications for strong effects at small quark masses, which have been predicted 
by ChEFT calculations. However, results at even lower quark masses with higher statistics will be 
required in order to confirm these predictions. 
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